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Abstract. We present a theoretical analysis of the evaporative cooling of a magnetically guided atomic 
beam by means of discrete radio-frequency antennas. First we derive the changes in flux and temperature, 
as well as in collision rate and phase-space density, for a single evaporation step. Next we show how the 
occurrence of collisions during the propagation between two successive antennas can be probed. Finally, 
we discuss the optimization of the evaporation ramp with several antennas to reach quantum degeneracy. 
We estimate the number of antennas required to increase the phase-space density by several orders of 
magnitude. We find that at least 30 antennas are needed to gain a factor 10* in phase-space density. 

PACS. 32.80.Pj, 39.25.+k, 05.30.Jp 



1 Introduction 

Evaporative cooling ^ is so far the only technique al- 
lowing the achievement of quantum degeneracy in dilute 
gases. Several models of evaporative cooling of a trapped 
cloud of atoms have been studied. They are based either 
on approximating the cooling process as a series of trun- 
cations of the distribution function followed by rethermal- 
ization [2] , or on an approximate analytical solution of the 
Boltzmann equation 21- 

With the achievement of Bose-Einstein condensation 
(BEG), the possibility of realizing sources of coherent mat- 
ter waves has arisen. Those "atom lasers" 0] open the way 
to fascinating applications in atom optics and interferome- 
try. To date, all atom lasers have been achieved in a pulsed 
mode: coherent streams of atoms were extracted from a 
Bose-Einstein condensate until it was completely depleted. 
As a first step towards the realization of a cw atom laser, a 
continuous source of Bose-Einstein condensed atoms was 
created by periodically replenishing a condensate held in 
an optical dipole trap with new condensates 0. This ap- 
proach, in combination with an appropriate outcoupler, 
would lead to a cw atom laser. 

Alternatively, the authors of [H] study the evaporation 
of an atomic beam propagating in a magnetic guide. They 
use an approximate solution of the Boltzmann equation 
based on a truncated gaussian ansatz P]. This approach 
permits the establishment of a set of hydrodynamic-like 
equations which is solved numerically. This treatment how- 
ever can be made rigorous only for a one-dimensional evap- 
oration. It would describe correctly evaporation on a di- 
electric surface on which atoms can be adsorbed 
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From an experimental point of view, the propagation of 
a single packet of atoms has been realized in both macro- 
scopic and microscopic atom guides |51I^1(JII11II12II13II14I 

[TB] . More recently, a continuous magnetically guided beam 
was achieved, by directly injecting a beam of cold atoms 
generated by a moving molasses magneto-optical trap |16) 
as well as by feeding a magnetic guide periodically at a 
high repetition rate ^ ^^^^ step of evaporation on 
this continuous guided beam has even been carried out 
by means of a single radio-frequency antenna. A natural 
way to attain degeneracy consists in using several radio- 
frequency (RF) antennas with decreasing frequencies. It 
would be the analog in space of the time-dependent RF 
ramp used in standard BEG experiments. The major dif- 
ference lies in the fact that the evaporative cooling is 
ensured by successive cycles of evaporation followed by 
rethermalization. 



In this work, we present the corresponding discrete 
model of evaporation by successive antennas. The two- 
dimensional character of the confinement allows us to de- 
rive analytical expressions that can be directly compared 
to experiments. After describing the magnetic guide in 
which atoms propagate, we derive how a single evapora- 
tion step affects the parameters of the beam (flux, tem- 
perature) and we deduce how the elastic collision rate 
within the beam, as well as the on-axis phase-space den- 
sity, evolve. We then show how one can characterize the 
occurrence of rethermalization between two successive an- 
tennas as was demonstrated experimentally in |17j . We 
finally address the issue of the optimization of the evapo- 
ration ramp with many antennas. 
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2 Magnetic guide 



3 Evaporation with one antenna 



Quadrupole guides can be produced by means of four wires 
or tubes equally spaced on a cylinder of radius a and car- 
rying currents +1 and — / alternatively. The cylindrical 
symmetry axis is chosen to be the z-axis in the following. 
The resulting magnetic field is well approximated by a lin- 
ear form: B{r) = br where r = {x^ + y^)^/^ is the distance 
from the z-axis and 6 = 2/.Jo//(7ra^) is the transverse gra- 
dient. In this quadrupolar configuration, atoms with low 
angular momentum are not stable against spin flips. To 
counteract this loss mechanism, one superimposes a bias 
field Bq along the axis of the guide ^Hl- As a consequence, 
the potential experienced by the low-field seeking atoms 
is: 

C/g(x,2/,z)=//[i?2 + 5V]i/2, (1) 

where /i is the magnetic moment of an individual atom. 

In the following, we consider a beam of atoms of mass 
m, transversally confined by the potential with a mean 
velocity -D, a temperature T, and a flux <P. 

At equilibrium the thermal average of Ug is: 



{Ug)T 



1 + a 



(2) 



In order to perform evaporation on the beam, one uses 
an RF antenna to produce an oscillating magnetic field 
at a frequency j^rf. Atoms whose trajectories cross the 
cylinder of axis z and radius R, defined by /ij^rf = Ug{R), 
will be resonant with the RF field and undergo a spin 
flip transition to an untrapped state. They are therefore 
removed from the atomic beam, whose flux decreases to 
the new value <P' . 

We define the evaporation parameter rj as rj = /i(^'rf — 
vo)/k^T, where vq = Ug{0)/h is the frequency correspond- 
ing to the bottom of the trap. In the following wc first 
calculate, as a function of rj, the fraction 0(7/) = <P' /ip 
of remaining atoms after such an evaporation cycle. The 
distribution function of the atoms just after evaporation 
is out of equilibrium. Assuming subsequent rethermaliza- 
tion of the beam, we then calculate its new temperature 
T' . In this paper we assume that the collision rate -Tcoii is 
much smaller than the typical period of oscillation (radial 
coUisionless regime). In other words, we assume that par- 
ticles do not undergo collisions over the range of efficiency 
of a radio frequency antenna. Therefore, determining if 
a particle will be evaporated or not depends only on the 
characteristics of its trajectory. 



where (3 ~ and a = ^.B^/UbT . For rubidium 87 

polarized in the low-field seeking state of the lowest hy- 
perfine level \F = —mp = 1), ^ = /j,_b/2 where jiB is the 
Bohr magneton, and a = 1 with _Bo = 1 Gauss for a tem- 
perature of T = 34 /xK. 

In the context of evaporative cooling, two quantities 
are of interest: The on-axis phase-space density V and the 
elastic collision rate I^-oii. For the potential we readily 
obtain 



1 1 <P f nb 



27rl + aw \kBT J {iirmkBTf/^ 



(3) 



Assuming that the clastic collision cross-section a is con- 
stant the collision rate reads: 



BcoW — 



1 + 2a ^' 



27r3/2 [l + af 



fib 
k^ 



knT 



(4) 



Two limits are of interest depending on the value of a. 
In the low temperature regime (a 3> 1), one can perform 
a Taylor expansion of resulting in the harmonic form 
with an energy offset [iBq: 



Ug{x,y,z) ~ fiBo + imcj^r^ 



(5) 



and angular frequency lu = [fj.b'^ /{mBo)]^^^ . The thermal 
potential energy (j^J reduces to fiBo + /cbT, and V (rcsp. 
Tcou) scales with # and T as <PT~^/'^ (resp. <PT-^/'^). In 
the opposite limit, where a <^ 1 one deals with a linear 
potential: Ug{x,y, z) c± iibr. The thermal potential energy 
is {Ug)j, = 2A:bT, and the scalings for V and r'coii are 
#T-7/2 and <PT-^/^. 



3.1 Fraction 0(77) = ^'/^ of remaining atoms 

To calculate ^(r/), we need to know the fraction of atoms 
in the beam whose trajectories in the [xy) plane cross the 
circle of radius R. 



3.1.1 Effective potential 

The confining potential Ug being cylindrically symmetric, 
the z component L of angular momentum is conserved. 
To study the radial motion of an atom, one therefore uses 
the effective onc-dimcnsional potential 



f/eff(L, r) = Ug{r) + 



2mr^ 



(6) 



which includes the centrifugal term. The radial coordi- 
nate r{t) of the particle oscillates between the turning 
points rmin and rmax- The atom is evaporated if and only 



if 1 



i? ^ r„ 



or equivalently if its transverse energy 



E fulfills E ^ Uos{R) (sec Fig.QJ. As a consequence, the 
evaporation criterion depends on the two constants which 
define completely the characteristics of the transverse tra- 
jectory: the total mechanical transverse energy E and the 
angular momentum L along the z-axis. 



3.1.2 Joint probability distribution of energy and angular 
momentum 

We now determine, for a thermal distribution at temper- 
ature T in the potential Ug, the joint probability distri- 
bution p{E, L) for an atom to have energy E and angular 
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Fig. 1. Effective potential Ucsir)- An atom with E ^ Ucb{R) 
is evaporated. 

momentum L. By definition 

p{E, L)^C j c-^'^ 5{E - H) 5{L - 4) dVd^p (7) 

where H(r,p) = /{2m) + Ug{r), = {rxp)-Cz and C is 
a normahzation constant determined by / p{E, L)dEdL = 
1. The joint probability can be readily recast in terms of 
the period of the radial motion Tr{E, L) [Hj: 

p{E, L) = C'c-P^Tr{E, L) . (8) 
3.1.3 Harmonic confinement 

For a harmonic confinement Ug{r) = ma;^r^/2, Tr{E, L) — 
tt/lu and E ^ En-ihi{L) = \L\uj. The latter condition just 
reflects that a finite angular momentum requires a min- 
imum energy which corresponds to the circular motion. 
We readily obtain the exact expression for the joint prob- 
ability: 

Ph.r{E, L) = ^^^^c-^^0{E - En,,n{L)), (9) 

where is the Heaviside step function. The fraction of 
remaining atoms is given by j2(J| : 

<^har(ry) = / Ph{E, L)dEAL = 1 - {^j^f'^c-^ . (10) 
Jt>(r) 

The integration domain in Eq. H1U|) is defined by 'D{R) = 
{Emin{L) = \L\u! ^ E ^ Ues{R)} and corresponds to the 
shaded area in the plane {E, L) depicted in Fig. |21 The 
subregion of r>(i?) such that L ^ Li{R) (resp. L ^ Li{R)) 
corresponds to trajectories not evaporated since r^ax < 
R (resp. Tmin > i?). The sohd line in Fig. 13 represents 
0hai (??)- For small i? (?7 ^ 1), few atoms have a sufficiently 
low angular momentum to be evaporated and ^har 1- 
For large R {rj ^ 1) the evaporative loss is negligible be- 
cause of the exponential decay of the energy distribution 
and (/)har ^ 1- Between those two limits, the fraction of 
remaining atoms has a minimum for rj = Tyj^^^j. = 0.5 with 

0har«ar) " 0-24. 



(a) E U.aiR) 




L 



-Li{R) Li{R) 
(b) E 




-Li{R) Li{R) 



Fig. 2. Domain of integration 'V{R) in the plane (E, L). Li{R) 
is the angular momentum of a particle having a circular tra- 
jectory of radius R. (a): harmonic confinement (with Li{R) = 
mR?uj), (b): linear confinement (with Li{R) — \/m^hR^). 

3.1.4 Linear confinement 

For linear confinement, the radial period Tr{E, L) has 
no simple analytical expression. However, Tr{E, L) has 
a weak dependence on the angular momentum L: for a 
given E^ when L varies from (linear trajectory) to its 
maximal value (circular trajectory), varies by less than 
10%. We therefore make the approximation Tr{E, L) ~ 
Tr{E,0) oc i?^/^ and approximate the joint distribution 
p{E, L) by: 

pUE. L) ~ CE^/^c-P''0{E - E^M) . (11) 

The minimum value of E for a given L (energy of circu- 
lar motion) is now: E^in{L) = 3{fib\L\)'^^^ / {2m^/^) . The 
fraction of remaining atoms can be written as (jjuniv) = 
Jx>(ii)Pl{E, L)dEdL where the integration domain is of 

the form V{R) = {Smin(i) ^ < Uos{R)} (see Fig. 

The result of the numerical integration based on the 
approximation (|11() is plotted as a dashed line in Fig. 
13 The shape is qualitatively the same as that for the 
harmonic confinement, but the minimum now occurs at 
V = '7rin - 1-25 with (^iin(r?i*„) ~ 0.33. 

To check the validity of the approximation (|ll(l . we 
performed a Monte-Carlo sampling of the atomic distri- 
bution in a linear potential, followed by elimination of 
particles whose trajectories crossed the circle of radius R. 
The result for 0iin(?7) is shown by dots in Fig. O The 
agreement is excellent. 
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Fig. 3. Fraction 4i{ri) = $'/<!> of remaining atoms after evapo- 
ration by one antenna, as a function of the evaporation param- 
eter rj. Solid line: Harmonic confinement; dashed line: Linear 
confinement, </!> being calculated with the approximation Ulljl : 
dots: Linear confinement, (fi being calculated by a Monte-Carlo 
simulation. 

It is convenient in practice to have a simple form for 
the expression of (punir]). Following the functional form of 
lfTU|l . we fit the curve </)iin('7) by a function of the form 
1 — arf'e~'^'^ where a, h, c are adjustable parameters and 
find a = 1.65, b = 1.13 and c = 0.92. The curve obtained 
this way differs from 0iin(?y) by less than 2%. This simple 
expression is used in |17| to fit the experimental data, with 
the temperature T of the beam being the only adjustable 
parameter. 

3.2 Temperature change T'/T 

Consider N atoms undergoing evaporation by one an- 
tenna. After evaporation the total transverse energy of 
the N' = (j}{ni)N remaining atoms is given by: 

{E±)^N f Ep{E, L) AEdL = ^Nk^Tgir^) , (12) 

J-D(R) 

where ^ = 2 for a harmonic potential and ^ = 3 for a linear 
potential. The dimcnsionlcss function g(r/) introduced in 
Eq. I|12|l can be calculated analytically for the harmonic 
confinement: 

5har(r/) = l-^^(3 + 2,7)e-''. (13) 

For a linear confinement, gwni'n) is calculated numerically 
thanks to the approximation Hll|l . 

The total (transverse and longitudinal) energy of the 
remaining atoms is E' = {Ej_) + N'kBT/2 + N'mv'^/2. 
After rethermalization (which occurs at constant total en- 
ergy), one has an equilibrium state with a temperature 
T' and energy E' = + l/2)N'kBT' + N'mv^/2. The 
relative change of temperature during this cvaporation- 
rethermalization cycle is therefore: 

r(^)---^fl + 2^44V (14) 
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Fig. 4. Temperature change r = T'/T after evaporation by 
one antenna and rethermalization, as a function of 77. Solid line: 
harmonic confinement; dashed line: linear confinement. 

In Fig. ^ T{ri) is plotted as a function of 77 (solid line: 
harmonic case, dashed line: linear case). For low 77 (ap- 
proximately lower than the value 77* for which reaches a 
minimum), T' > T: since essentially low energy atoms are 
evaporated, the remaining beam acquires a higher temper- 
ature after thermalization. The opposite occurs for high 77: 
the beam is cooled by evaporation of high energy particles. 



3.3 Variation of collision rate and phase-space density 

After an evaporation followed by a subsequent rethermal- 
ization, the flux is multiplied by a factor 0, and the beam 
temperature is T' . Wc can then calculate, using Eqs. © 
and the new phase-space density V and the new col- 
lision rate -TpQn. 

The variation in phase-space density ^(77) = T>' /T> is 
plotted as a function of 77 in Fig. [3 (solid line: harmonic 
case; dashed line: linear case). For a linear confinement an 
increase of V is obtained as soon as 77 > 2, and for rj ~ 3.0 
a maximum gain in phase-space (by a factor 1.86) can be 
achieved with a single antenna. The gain is less significant 
in the case of a harmonic confinement, reflecting the less 
favorable scaling law. 

In order for evaporative cooling to be efficient, the elas- 
tic collision rate must not decrease during the process, 
otherwise the time needed for rethermalization becomes 
prohibitively long. Fig. depicts the variation of colli- 
sion rate 7(77) = /^coii/^coU for one evaporation cycle. The 
striking result is that in the case of harmonic confinement 
the collision rate cannot increase significantly (7 reaches a 
maximum value of 1.0015 for 77 = 6.5). This behaviour is 
due to the two-dimensional geometry of the problem and 
does not hold for a three-dimensional harmonic confine- 
ment. In contrast, for a linear guide, and for 77 sufficiently 
high, the collision rate can increase by several percent at 
each cycle, creating the possibility of runaway evaporation. 
This salient feature of the linear potential is well known 
for three-dimensional evaporative cooling 
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Fig. 5. Gain in phase-space density D' I'D after evaporation 
by one antenna and rethermalization, as a function of r]. Solid 
line: harmonic confinement; dashed line: linear confinement. 
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Fig. 6. Variation of the collision rate after evaporation by one 
antenna and rethermalization, as a function of t). Solid line: 
harmonic confinement; dashed line: linear confinement. 



3.4 Multi-radii evaporation 

The previous evaporation method has the disadvantage 
that for any finite ?/, some atoms with large energy and 
angular momentum are not evaporated. The contribution 
of these atoms to the energy of the truncated distribution 
is not negligible, and therefore the cooling efficiency of the 
evaporation is limited. 

In order to improve this efficiency, one can evaporate 
all the atoms whose trajectories lie, at least partially, out- 
side of the cylinder of radius R. To realize such an evap- 
oration, one must use a "continuum" of evaporation radii 
Re G c)o]. In practice, this can be achieved by scanning 
the radio- frequency from Ug{R)lh to a value large 
compared to ksT/h, above which the population of atoms 
is exponentially small. The scanning rate must be large 
compared to the inverse of the time t^v spent by an atom 
in the range of efficiency of an antenna, and small com- 
pared to the inverse of the mean radial period Tr{E, L). 
For the typical parameters of the experiment described 



Fig. 7. Fraction = <!>' /<1> of remaining atoms after evap- 
oration by one antenna, as a function of 77, for the multi-radii 
evaporation scheme. Solid line: harmonic confinement; dashed 
line: linear confinement. 

in |17j . a scanning frequency of 50 Hz would fulfill both 
criteria. 

One can calculate, in a way similar to the calcula- 
tions of the previous subsections, the change in atomic 
ffux (f){r]) = in temperature riji) = T' /T, in phase- 

space density 15(77) = T)' /T), and in collision rate 7(77) = 
^cou/^coiii ^^^^ improved evaporation scheme. The do- 
main of integration is that shown in Fig. [21 but reduced 
to the angular momentum lower than Li{R). In the case 
of harmonic confinement, the calculations can be done ex- 
plicitly and yield the following results: 

0(77) = 1 - e-^-) _ (7r77)i/2e-"erf (77^/2)^ 
.9(77)= 1 -(l + 77/2)e-2" 

- (3/4 + 7;/2)(7r77)i/2e-"erf (77^/2) (^5) 

where erf(x) = 2tt^^^^ exp{—t^)dt is the error func- 
tion. The complete results for 0, r, 7 and S and for both 
types of confinement are plotted in Figs. [3 to EH Obvi- 
ously, with this new evaporation scheme, the number of 
atoms remaining vanishes when 77 — > 0. The temperature 
change (Fig. IHJ is now a monotonous function of 77 and for 
77 — > 0, in the harmonic (resp. linear) case one has r — > 1/5 
(resp. 1/7) since only the longitudinal kinetic energy con- 
tributes to the energy of the remaining particles. 

For a linear potential, the gain in phase-space density 
without decrease in collision rate is 2.56 (at 77 = 2.5) as 
compared to 1.86 (at 77 = 3.0) in the single radius evapo- 
ration scheme. 



4 Evaporation with two antennas 

The presence of elastic collisions in the atomic beam can 
be probed with a two-antenna experiment |17| . A first an- 
tenna with frequency vi evaporates the atoms with an 
evaporation parameter 771 = hvi/ [k^T), where T is the 
initial temperature of the atoms. The remaining ffux is 
<!>' = (f>{r]i)<P. A second antenna with frequency 1^2 ~ 
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Fig. 8. Temperature change r — T' /T after evaporation by 
one antenna and rethermalization, as a function of rj, for the 
multi-radii evaporation scheme. Solid line: harmonic confine- 
ment; dashed line: linear confinement. 




Fig. 9. Gain in phase space density after evaporation by 
one antenna and rethermalization, as a function of rj, for the 
multi-radii evaporation scheme. Solid line:harmonic confine- 
ment; dashed line: linear confinement. 



1.25 




Fig. 10. Variation of the collision rate after evaporation by 
one antenna and rethermalization, as a function of rj, for the 
multi-radii evaporation scheme. Solid line:harmonic confine- 
ment; dashed line: linear confinement. 



r]2kBT/h, placed at a distance d downstream, is used to 
probe the distribution of atoms after this first evaporation 
stage. The flux after this second antenna is denoted <P" . 
In the absence of collisions between the two antennas, the 
distribution remains out of equilibrium. In particular, if 
1^2 ^ 1^1, the second antenna cannot remove extra atoms, 
since all atoms whose energy and angular momentum cor- 
responded to the evaporation criterion have already been 
removed by the first antenna. But even a partial rether- 
malization leads to some extra losses from the second an- 
tenna, since clastic collisions provide a redistribution of 
energy and angular momentum. If the beam completely 
rethermalizes during its propagation between the two an- 
tennas, we find a spectrum similar to that of Fig. |21 reveal- 
ing the new temperature T' . 



4.1 In absence of collisions 

We calculate here the fraction <P" /<P' in the limit where 
rethermalization does not occur between the two anten- 
nas, i.e. when the inequality r'coii-'^ <C w is fulfilled. After 
evaporation by the first antenna at rji (corresponding to 
an evaporation radius Ri), the out of equilibrium joint 
probability distribution pi{E, L) of {E, L) is given by 



PiiE, L) = p{E, L) X e([/eff(i?i) - E) . 



(16) 



It is non-zero only in the domain Vi U X'2 of Fig. ^2 The 
fraction of remaining atoms <^"/<?' is therefore given by 
integrating over the domain P2 = {Enun{L) ^ E ^ 
min(C/eff(i?i),t/cff(i?2))}: 



„ [ pi{E,L)dEdL 



(17) 



pi{E, L)dEdL 
V1UV2 



The result of the numerical integration is shown in Fig. ^] 
as a dashed line, for the specific case of a linear confine- 
ment with rji = 2. As expected, the fraction <P" /<P' reaches 
1 for r]2 = 2. 



4.2 Rethermalization 

If rethermalization occurs over a distance d, that is if the 
collision rate is such that: r'coii-c' ^ v, the second antenna 
evaporates a beam of temperature T' = T{rii)T at thermal 
equilibrium. The fraction #"/<P' is therefore: 



hi>2 



V2 

t(t?i) 



(18) 



This curve is plotted as a solid line on Fig. ^| The large 
difference between the collisionless and collisional regimes 
makes it easy to detect experimentally, by scanning V2 in 
the vicinity of vi , the presence or the absence of collisions 
within the beam )17| . 
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(a) E/k^T 
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(b) E/kBT 



L 



Fig. 11. The integration domains in the plane {E, L) involved 
in the calculation of in the collisionless regime, (a): Har- 

monic confinement; (b): Linear confinement. 
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Fig. 12. Fraction <1>" of the flux remaining after the second 
antenna as a function of its parameter rj2- The first antenna 
is driven at a frequency v\ — rjiksT/h with rji = 2 = I.677*. 
The dashed line corresponds to the collisionless regime; the 
dotted line to a mean number of collisions between antennas 
Nc = 1.1; the solid line to full rethermalization. These curves 
are calculated with a numerical simulation, in the case of a 
linear confinement. 



4.3 Number of collisions required to thermalize 

In order to give a quantitative theoretical account for the 
distance d needed to rethermalize in a two-antenna exper- 
iment, we have developed a numerical simulation of the 
motion of the atoms in the magnetic guide for both linear 
and harmonic confinements. We use a molecular dynamics 
simulation |23 taking into account the evaporation crite- 
rion. 




Fig. 13. Distance from equilibrium S{rj\,Nc) (see text) as 
a function of the mean number of collisions Nc between the 
two antennas, for different values of 771. These data have been 
obtained by a Monte-Carlo simulation performed for a linear 
confinement. The relaxation towards equilibrium is clearly ex- 
ponential for r]\/rf = 0.3, and approximately exponential for 
771/77* = 3.2. For 771/77* = 1.6, the first antenna puts the gas in 
a strongly out-of-equilibrium state, and the relaxation is not 
exponential. 



The simulation is performed in dimcnsionlcss units. 
The unit length r„ is defined by Ug{ru)~Ug{Q) = {{Ug)T-' 
J7g(0))/2, the distance between the two antennas is d = 
5000 r„. The mean velocity is 10 times the thermal veloc- 
ity (fesT/m)^/^. The phase-space variables of each particle 
are evolved by advancing the position and the velocity at 
discrete time steps At according to a second-order sym- 
plectic integration |22| . Using the cylindrical symmetry we 
restrict the evolution to the first quarter (a; > 0, 2/ > 0) 
of space with reflecting walls at planes x = and y = i). 
Binary elastic collisions are taken into account using a 
boxing technique The constant cross section used to 
calculate the probability of collisions is adjusted to the 
desired number of collisions between the two antennas. 
Simulations have been performed with 1.5 x 10^ particles 
and 5.4 x 10^ boxes. The time step At is chosen to be 
small with respect to the typical collision timcscale and 
the typical period of oscillation. 

We simulate a two-antenna experiment in which the 
first antenna is operated at a fixed 771 . The mean number 
of collisions per atom during the propagation between the 
two antennas is denoted N^. The second antenna is op- 
erated at 772, and we calculate <?"/<?' as a function of 772. 
An example of such a curve, for a linear confinement, and 
with Nc = 1.1, is plotted as a dotted line in Fig. 1121 

We then infer the ratio f{rii,Nc) = <?"/<?' (772 = 771) 
as a function of Nc- For a collisionless propagation {Nc = 
0, dashed line in Fig. I12() . wc have f{r]i,0) = 1. For a 
complete rethermalization, {Nc ^ 1, solid line in Fig. 112(1 . 
the ratio tends to the limit foo{'>li) = 

For a linear confinement, we depict in Fig. 1131 the nor- 
malized quantity 

^( ,r N f{Vi,Nc) - /oo(77i) 



8 



T. Lahaye and D. Guery-Odelin: Discrete-step evaporation of an atomic beam 



which measures the distance from equihbrium, as a func- 
tion of the mean number of colhsions Nc- For a given rji 
and in the absence of coUisions Nc = 0, S = 1. The re- 
laxation of S towards zero with increasing number of col- 
hsions Nc reflects the rethermalization process. For 771 < 
0.377* or 771 > 2.577*, S decays approximately in an ex- 
ponential manner with Nc'. ^ — exp[ — A^c/^c(^i)]- 
example we find 77c(4) = 2.5 ± 0.5 for the linear confine- 
ment and 77c (2) = 1.3±0.2 for the harmonic confinement. 
For intermediate values of the ratio 771/77*, we obtain a 
non-exponential decay (see Fig. E| for 771/77* = 1.6). In 
this range 0.3 < 771/77* < 2.5, more than 50% of the atoms 
are evaporated by the first antenna, leading to a state of 
gas that is very far from equilibrium. 

This calibration by numerical simulation allows one to 
determine the collision rate by following the thermaliza- 
tion process for different distances between the two anten- 
nas [Tt] . 



5 Evaporation with many antennas 

The realization of a magnetically guided atomic beam in 
the collisional regime has recently been performed |17) . 
The phase-space density of this beam is on the order of 
2 X IQ-®. To gain the 8 orders of magnitude required to 
reach degeneracy one can implement evaporative cooling 
by means of several successive antennas. In this section wc 
assume that the frequency of each antenna is such that the 
parameter i] is constant throughout the evaporation. 

The number of antennas Af{r]) required to achieve a 
gain A in phase-space is, for a given value of 77, M(ri) = 
ln(Z\)/ ln(i5(77)). We want to minimize the guide length 
needed to reach degeneracy. This length £(77) can be eval- 
uated from the initial collision rate of the beam r^c,n- 
twcen the antennas n and n+1, the collision rate is r^^^^j". 
If the number of collisions required for full thermalization 
between two successive antennas is A^o (^0 ^ "-c('7)), the 
distance between antennas n and 77+1 has to be larger 
than vNo/{^°^^^Y')■ Therefore the total length is 



vNq 

^ coll 



1 + 7 + ... +7' 



-A/'+l 



vNo 1 - 7-^ 



r°„ 1 



T 



(19) 



For a given value of A, there exists an optimal choice 
77 for the parameter 77 such that £ is minimized. For this 
optimization we assume that A^o does not depend on 77. 
The results are shown, for both types of confinement and 
for the two evaporation schemes, in Tabled The evapo- 
ration is much less efficient for the harmonic confinement, 
which reflects the absence of runaway evaporation for this 
type of trap. For the linear confinement, one sees that the 
multi-radii evaporation scheme leads to a slightly shorter 
evaporation length, at the expense of a reduced output 
flux. However a more careful study shows that for com- 
parable output fluxes, the two schemes require approxi- 
mately the same length and antenna numbers. In the case 



Scheme 


Single radius 


Multi-radii 


Trap 


Har 


Lin 


Har 


Lin 




211 


12.0 


211 


10.5 


V 


4.6 


4.3 


4.6 


3.6 


N 


142 


42 


142 


29 




3.9.10-3 


5.4.10"* 


3.9.10"3 


1.3. 10-* 


Tf/T, 


6.9. 10"^ 


6.0.10"* 


6.9. 10"^ 


3.8. 10-* 



Table 1. Parameters of the optimized (minimum length) evap- 
oration ramp for harmonic (Har) and linear (Lin) confine- 
ments, for the two schemes of evaporation, with a total gain 
in phase-space density of Z\ = 10*: dmin, evaporation length 
(in units of vNq/ rj.'^n); fj, parameter 77 for minimal evapora- 
tion length; N, number of antennas needed; ^f/'l'i, total fiux 
reduction due to evaporation; Tf/Ti, total temperature reduc- 
tion due to evaporation. 



of harmonic conflnement, since the collision rate cannot 
increase significantly, minimizing the evaporation length 
requires that one operates at a large 77 (where 7 ~ 1). For 
such high 77, the difference between the two evaporation 
schemes is negligible, as can be seen in Table Q] 

In this paper, two competing effects that should be 
considered in a more realistic model of the evaporation 
ramp have been neglected. First, the real shape of the 
semi-linear potential ^ should be taken into account. In- 
deed, even if at the beginning of the evaporation ramp 
the temperature is high enough for the linear potential 
to be a good approximation (a <C 1), as the temperature 
is reduced the effect of the bias field cannot be neglected 
anymore when a > 1. This effect may lead to a longer 
evaporation length since, for the harmonic confinement, 
the evaporative cooling is less efficient. The second effect, 
neglected by the classical description used in this paper, is 
the bosonic stimulation close to degeneracy, which might 
enhance the efficiency of the last evaporation steps. 



6 Conclusion 

We have studied theoretically the discrete-step evapora- 
tive cooling of a magnetically guided atomic beam. First, 
the action of a single antenna on the beam is analyzed in 
detail. Then we clarify how two-antenna experiments can 
lead (i) to the characterization of the collisional regime 
of the beam and (ii) to the quantitative estimate of the 
collision rate, by studying the number of collisions needed 
for thermalization of the truncated distribution function. 
Finally, the problem of optimizing the evaporation ramp 
with many antennas is addressed. We find that gains in 
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phase-space density of 10^ require a total of at least ^ 30 
evaporation steps. 
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P. Cren for the development of a preliminary version of the 
numerical simulation. We acknowledge stimulating discussions 
with J. M. Vogels, K. J. Giinter, K. Kozlowski and the ENS 
laser cooling group. We thank D. Roberts for careful reading 
of the manuscript. This work was partially supported by the 
CNRS, the Region Ile-de- France, the French Ministry of Re- 
search, the College de France, the Ecole Normale Superieure, 
the Bureau National de la Metrologie, the University of Paris 
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